Genuine tripartite entanglement in the non-interacting Fermi gas 
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We study genuine tripartite entanglement shared among the spins of three localized fermions in 
the non-interacting Fermi gas at zero temperature. Firstly, we prove analytically with the aid of 
entanglement witnesses that in a particular configuration the three fermions are genuinely tripartite 
entangled. Then various three-fermion configurations are investigated in order to quantify and 
calculate numerically the amount of genuine tripartite entanglement present in the system. Further 
we give a lower and an upper limit to the maximum diameter of the three-fermion configuration 
below which genuine tripartite entanglement exists and find that this distance is comparable with 
the maximum separation between two entangled fermions. The upper and lower limit turn to be 
very close to each other indicating that the applied witness operator is well suited to reveal genuine 
tripartite entanglement in the collection of non-interacting fermions. 

PACS numbers: 03.65.Ud, 03.67.Mn, 71.10.Ca 



I. INTRODUCTION 

Entanglement is in the heart of quantum mechanics 
and of great importance in quantum information theory. 
Two entangled particles already offer a valuable resource 
to perform several practical tasks such as quantum tele- 
portation, quantum cryptography or quantum computa- 
tion However, the multipartite setting due to the 
much richer structure suggests many new possibilities 
and phenomena over the bipartite case. Indeed, mul- 
tipartite states may contradict local realistic models in 
a qualitatively different and stronger way 0] • Moreover, 
this feature allows to implement novel quantum infor- 
mation processing tasks such as quantum computation 
based on cluster states entanglement enhanced mea- 
surements [4|, quantum communication without a com- 
mon reference frame [5| and open-destination teleporta- 
tion @. 

Though the characterization of multipartite entangle- 
ment is studied in great depth @> H|i still rarely in- 
vestigated in solid state systems (see 9| and references 
therein), however, it definit ely play s an essential role in 
quantum phase transitions [10 . [ll| and might well be a 
key ingredient to unresolved problems in physics such as 
high temperature superconductivity [T2"j . In this article 
we investigate genuine multipartite entanglement shared 
among the spins of three fermions in the Fermi gas of 
non-interacting particles at zero temperature (degener- 
ate Fermi gas) following the work of Refs. [H, Il4 | We 
apply an entanglement witness developed in Ref. [9( (and 
also appeared in Ref. [l5| ) in order to reveal genuine tri- 
partite quantum correlations in the collection of fermions 
and using the approach in Ref. [l6| we also characterize 
quantitatively the amount of it. 

The article is organized as follows: In Sec. ITI1 according 
to Refs. [HI, 03] we present the three-spin reduced den- 
sity matrix of the degenerate Fermi gas, and show by an- 
alytical means (extending the related results of Ref. fl4l |) 
that both GHZ-type and W-type witnesses are unable 



to detect genuine tripartite entanglement (GTE) among 
the spins of three localized fermions. On the other, we 
demonstrate that for specific configurations of the three 
fermions the GTE witness of Ref. [9|] is capable to signal 
GTE both for the two- and three-dimensional (2D and 
3D) degenerate Fermi gases. In Sec. Mil on the basis of 
this witness a formula is constructed to the lower bound 
of the generalized robustness (En) of genuine tripartite 
entanglement. With the aid of this formula in Sec. II VI we 
quantify numerically genuine tripartite quantum correla- 
tions for various arrangements of the three particles. In 
Sec. [V] we determine lower and upper bounds to the GTE 
distance (i.e., to the largest diameter of the three-fermion 
configuration below which GTE is still present in the sys- 
tem) both in the 2D and 3D degenerate Fermi gases. The 
paper concludes in Sec. I VII with a brief summary of the 
results obtained and discusses possible schemes to extract 
GTE from the system. 



II. ANALYSIS OF THE DENSITY MATRIX FOR 
THREE FERMIONS 

A. Three-spin reduced density matrix 

Consider a system of N non-interacting fermions in a 
box with volume V. At zero temperature the ground 



state of the system is \4>q) = U k F c], a \vac) , where fcp = 

(2>tt 2 N /V) 1 ^ 3 is the Fermi momentum and \vac) denotes 
the vacuum state. From the ground state of the system 
\cf>o) one obtains the three-spin reduced density matrix 
(up to normalization) between the three fermions local- 
ized at positions r, r', and r", 



P3(s, 



';t,t',t") = 



(0o|^(r)^(rO^(r'')^"(r'')^(r , )^(r)|0o) , (1) 

where V's(r), V4( r ) are ne ^ annihilation/creation opera- 
tors for a particle with spin s located at position r sat- 
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isfying {-0 s (r), -0j(r)} = 5 StS 'S(r — r'). The justification, 
that the above matrix elements indeed describe three- 
qubit quantum states is discussed in the Appendix of 
Ref [l||] . Following Refs. [H, 03 the explicit formula for 
the three-spin reduced density matrix p$ is given by 

pa = (i-p)^+pi2i*r 2 >(*r 2 i® y +pi3i*r 3 >(*r 3 i® ^ 



■P23|* 23 >(* 



2:1 
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(2) 



- P13 + P23 and I„ 
Further, <EC 



denotes the n x 



where p = P12 

identity matrix. Further, 4C. = (| |4> - I IT))/V2 is 
the singlet state of the pair ij in the orthonormal basis 
{| t), I 4)}. The value p^ depends only on the relative 
distance between the three fermions and can be written 
explicitly for the fermion pair ij as [l4j 



P, :! 



fij ~l~ fij fikfjk 



2 + fij + ffk + fjk fij fikfjk 



(3) 



where the analytic form of fij depends on the spatial 
dimension of the system, that is we may write 



flf = 2J 1 (k F r lj )/k F r l j 



(4) 



in the case of the two- and three-dimensional Fermi gases 
[l9l | . In the above formulae j\ and Ji denote the spher- 
ical and the first order Bessel function of the first kind, 
respectively. 

Actually, owing to the collective SU(2) rotational 
symmetry of the model Hamiltonian of non-interacting 
fermions, many matrix elements of p% in © are forced 
to be zero. Explicitly, the states which are invariant un- 
der collective SU (2) rotation of the three qubits are the 
three-qubit Werner states [2(J, and they can be given in 
the form [H 



P : 



ET_k_ 
4 

Ai= + ,0,1,2,3 



(5) 



where Rk are certain linear combinations of permutation 
operators and rk (p) — Tr(pRk)- Using the definitions for 
Rk from Ref. [2(| and the explicit form of the state p$ 
from 10), we are able to calculate the parameters rk for 
the three-spin reduced density matrix p 3 , which read as 
follows 

l-p 



n 



2 

1+p 
2 

P12 + P13 - 2p 



2:', 



Ol3 -Pl2) 



2V3 



r 3 = 



(6) 



B. Possible range of parameters pij 

According to the Lemma 2 of Ref. [2(| /J in (J5J) is a 
density matrix only if r+, ro > 0. These inequalities im- 
ply for the state P3 by the virtue of (J6j) that p lies in the 
interval 



1 < P < +1 



(7) 



Let us observe in ([3]) that \fij\ < 1 both for the 2D and 
3D Fermi gases. This fact together with the bound to p in 
([7|) and also the definition p = P12 + P13 +P23, after some 
algebraic manipulations (which are not detailed here), 
lead to the bounds 



1 < Pij < 1 



(8) 



for the three different fermion pairs ij — 12, 13, 23. Thus, 
the parameters p^ appearing in state ps are limited by 
the values ±1. 

Let us introduce the class of biseparable three-qubit 
states B, i.e., the states 



P : 



1 P l \^i){^i 



(9) 



which can be expressed as a convex sum of projectors 
onto product and bipartite entangled vectors [7[. In 
definition §§§ the pure states 1^) are separable on the 
Hilbert-space of three qubits 1,2,3 with respect to one 
of the three bipartitions 1|23, 12 1 3 or 13 1 2 and Pi > 
adding up to 1. We say that a general three-qubit state is 
genuine tripartite entangled when it is not in the class of 
biseparable states B, that is they cannot be constructed 
by mixing pure states containing bipartite entanglement 
at most. Clearly, if p^ in ((2]) was positive for each of the 
three different pairs (whose sum p is upper bounded by 
+1 according to (O), then p 3 should define a biseparable 
state. 

Note, that the bound to Pi/s in ([8|) may not be tight, 
therefore it is not evident whether py can take up neg- 
ative values at all. However, by arranging the three 
fermions in a particular geometry, we demonstrate that 
Pij may take the value — 1/3 as well: Taking the first and 
second derivatives of the functions fij(kpr) (defined for 
both the 2D and 3D Fermi gases under equations (|4|)) 
with respect to x = kpr, we observe that in the limit 
x — > they behave as 



lim f^ (a:) = 1, lim f[ 2 (x) = 0, 



hm/£Or)^0 

IE— >0 J 



(10) 



Now let us place the three particles on a line so that 
particle 2 would lie just at the midpoint between particle 
1 and particle 3, and let the relative distance between 
these outer particles tends to zero. Actually, in the limit 
x — > 0, Pij in © can be given explicitly by applying 
PHospital's rule twice and by taking account the limiting 
values (flOl). As a result we obtain 



P12 = P23 = 2/3, 



P13 



-1/3 



(11) 
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both for the 2D and 3D Fermi gases. 

In the next two subsections we propose witness oper- 
ators in order to reveal GTE in the degenerate Fermi 
systems. An observable which is well suited for signal- 
ing genuine tripartite quantum correlations in Heisenberg 
spin lattices |9[] turns out to detect GTE in the degener- 
ate Fermi gas as well. 



C. Generalized GHZ and W witnesess 

For deciding whether the state p 3 with the explicit pa- 
rameters P12 — P23 — 2/3, P13 = —1/3 in (fTTj) is genuine 
tripartite entangled, we will use entanglement witnesses. 
A witness of genuine tripartite entanglement is an ob- 
servable IT with a positive mean value on all biseparable 
states so a negative expectation value Tr(/TI) guarantees 
that the state p carries genuine tripartite entanglement 
[2ll ]. Thus a witness operator which separates genuine 
tripartite entangled states from the biseparable set B 
(defined by equation (J9j) ) can be given in the form [22j 



where 



W^ = Al8-|^| , 



A= max |(#A>| 2 

\<t>)€B 



(12) 



(13) 



A simple method has been found in Ref. [22| to deter- 
mine A for any pure genuine tripartite entangled state 
In particular, let be the GHZ-\ike and the W- 
like state [23[, which are respectively 



\GHZ(a)) 



|W(/?, 7 )) = 



|ni,n 2 ,n 3 ) +e M | -ni,-n 2 ,-n 3 ) 



Piji ij — 12, 13, 23, hence no GTE could be revealed in 
the three-spin reduced density matrix p 3 by the applica- 
tion of these witnesses. 

We confirm and extend this result by generalizing the 
GHZ and W witnesses to the form (fTB"]) with the cor- 
responding values a, f3, 7 € [0, 2tt] and with the arbi- 
trary local bases |n,) = cos(^/2)| f) + e^sin(^/2)| | 
), i = 1,2,3. Then the trace of p 3 II, where II denotes 
either WcHZ(a) or Ww(/3,j) an d 93 is the state ([2]), is 
a linear combination of the trigonometric functions co- 
sine/sine with arguments 6i,4>i, i = 1,2,3. Owing to 
convexity arguments Tr(p 3 II) can be extremal only if 
Pij 6 {+1, —1}, ij = 12, 13, 23 in the permitted range ([8]) 
and 6i,^i £ {0, tt/2, tt, 3tt/2, 2tt}, i = 1,2,3. More- 
over, we may fix three parameters, e.g., (j>i = $2 — 
and 9\ = 0, owing to the invariance of the state p 3 un- 
der the collective SU(2) rotation. Considering all the 
possible combinations of the remaining parameters py, 
ij = 12, 13, 23 and 6* 2 , 03, 4>3 from the set above, we found 
that Tr(p 3 II) > is always true for the witness operators 
II in (I15p . Therefore, genuine tripartite entanglement 
could not be witnessed by the general GHZ-type and 
W-type witnesses (|15l) . as well. 

Naturally, we may ask whether there exist other wit- 
nesses over the GHZ/W- types in (fT5]) which are bet- 
ter suited for detecting GTE in the state p 3 . Indeed, in 
Ref. [23j it has been shown that there are several genuine 
triparite entangled states which are not witnessed by the 
operators (fT5|) . Note, that nonlinear entanglement wit- 
nesses [24j may show improvement with respect to linear 
witnesses in the multipartite case as well. However, if 
we stick to linear functionals we are even able to reveal 
GTE in the state p 3 , as it will be discussed in the next 
subsection. 



ni,n 2 , -n 3 ) + e'Hni, -n 2 , n 3 ) + e* 7 | - ni,n 2 ,n 3 



V3 



>D. The witness operator of Giihne et al. 



(14) 



where {n.;,— n^} denotes an arbitrary local orthonormal 
basis in the Hilbert space of qubit i. Note, that for the 
original GHZ state the phase a = and for the original 
W state the phases j3 = 7 = 0, and the corresponding pa- 
rameters A appearing in the witness operator (|12[) are 1/2 
and 2/3, respectively Q- The states (TT41 . however, can 
be transferred to the original ones by local unitary oper- 
ations, which leave the parameters A unchanged, that is 
we have the witness operators 



Giihne et al. [9] showed that the internal energy is 
a good indicator of genuine tripartite entanglement in 
macroscopic spin systems. The idea was to write the 
internal energy in terms of the mean value of the observ- 
ables Wijk — S % ■ + $ ■ a k , where a 1 = [a l x , a l y , <j l z ) 
is the vector of Pauli spin operators associated with the 
qubit i, and the absolute value of (W^k) has been shown 
to be a witness itself, capable to detect GTE. Namely, it 
has been proven @ that if the inequality 



W, 



GHZ(a) 



2 S 



\(W ijk )\ > 1 + V5~ 3.236 



(16) 



\GHZ(a))(GHZ(a)\ 



W w{0 ^ = -I 8 -\W(0, 7 ))(W(fr)\ 



(15) 



for the GHZ- and M^-like states, respectively. 

Lunkes et al. [l4| applying these witness operators for 
the special case a = (3 = 7 = and |n.;) = | |), i = 1, 2, 3 
established that neither Tt(p3Wghz) nor Tr(p 3 Ww) can 
become negative in the permitted range of parameters 



holds, the qubits k are genuinely tripartite entangled. 

Let us use this inequality (fTB|) in order to reveal GTE 
in the degenerate Fermi gas among the spins of the 
three fermions i,j,k. Plugging the state p 3 in @ into 
the expectation {Wijk) = Tr(p 3 Wyfe) one has {Wijk) — 
3(py +pjk), which by substitution back into (|16p gives 
the condition 



\{W ijk )\ -=^\Pij+Pjk\ > l + VE 



(17) 
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for the existence of genuine tripartite entanglement in the 
degenerate Fermi gas. Next we discuss from the view- 
point of witnessed GTE by the mean of this condition 
two different three-fermion configurations: 

(a) Consider the case investigated before in Sec- 
tion HlBl that three particles lie evenly spaced on a line 
close to each other. Choosing ijk = 123 and recall- 
ing pi2 = P23 = 2/3 from (fTTj) . by the virtue of (fT7j) . 
|(Wi23)| = 4 > 1 + \/5, hence the three-fermion state pz 
corresponding to this arrangement of particles is genuine 
tripartite entangled. 

(b) In this case the particles are separated from each 
other by equal distances, i.e., the particles are put on 
the vertices of an equilateral triangle. Owing to three- 
fold symmetry of this configuration the state p 3 contains 
an equal mixture of maximally entangled states l^ 1- ), 
that is, all Pij, ij — 12, 13, 23 in © must have the same 
value. Further, considering the constraint \p\ = \p\2 + 
P13 + P23 1 < 1 in © and also owing to the left-hand 
side of (ED we have |(W 12 3>| = I (W331) | = |(Wi32>| = 
2\p\ < 2 implying that in this case all possible KlFijfc)! 
(with different permutations of ijk) are smaller than the 
bound 1 + V5- Consequently, no GTE can be revealed by 
the witness (fl~6| of Guhne et al., no matter how far the 
fermions are separated from each other. It is reasonable 
to think that there is indeed no GTE associated with this 
highly symmetrical configuration, as it has been argued 
in Ref. [14| by attributing it to the Pauli principle. In 
the next section we construct from the observable Wijk a 
witness operator Wijk, which has a maximum eigenvalue 
smaller than unity (Wijk < Is) an d with the aid of it 
a lower bound is given for the amount of GTE in the 
state /?3 quantified by an entanglement monotone, the 
generalized robustness Er. 



The generalized robustness Er as a GTE measure 
quantifies how robust the genuine tripartite entangled 
state p is under the influence of noise, and also has a ge- 
ometrical meaning measuring the distance of p from the 
biseparable set B [25| . According to Ref. [16[ En can be 
expressed in a Lagrange dual representation 



Er(p) — max{0, — min Tr(/oII)} 



(19) 



where the set M is given by the restriction II < Is for 
the GTE witnesses II. 

Since Wijk with any permutation of ijk defines a valid 
GTE witness with maximum eigenvalue smaller than 
unity, we are able to develop the lower bound 

£fl,mta(p) =max{0,- min Tr{pW lok )} (20) 

ijfcG{123,231,132} 

to the generalized robustness (TiT)|) of an arbitrary state 
p on qubits 123, as it is discussed in Refs. [Ill [27]] ■ This 
lower bound by applying (fT8|) for the particular state ps 
in ^ reads as 



-ER,min(p3) = max {0, 

r/fce{123,231,132} 



3(pij +Pjk) -i-VE 

5 + V5 



}• 

(21) 

In the next section this formula will be applied to give 
explicitly a lower bound to the generalized robustness 
En of GTE for various configurations of three fermions, 
associated with the reduced state p 3 of the degenerate 
Fermi gas. 



IV. NUMERICAL CALCULATIONS TO THE 
LOWER BOUND OF E R 



III. DERIVING A LOWER BOUND TO THE 
GENERALIZED ROBUSTNESS E R 

Up to this point the observable Wijk was applied 
for witnessing genuine tripartite entanglement. On the 
other, they are also good for quantifying it [l6| (see [l7| 
as an application to a magnetic material). The maxi- 
mum eigenvalue of the operator —Wijk is 4, thus consid- 
ering ([TB)1 we may construct the following witness oper- 
ator, 



Wijk = 



(i + Vsjh - w tJk 

5 + V5 



(18) 



whose negative mean value Ti(pWijk) guarantees that 
the three-qubit state p is genuine tripartite entangled. 
Further the witness operator is normalized so that 
Wijk < Is- R is apparent that provided the mean value 
(Wijk) > for state pz (i.e., p^ + Pj k > according 
to calculations in Sec. IIID[) . the witness operator (Tj"gj) is 
just as powerful to detect GTE associated with state p$ 
as KlFijfe)! in the inequality (QI 



A. Fermion moving on a straight line 

Now we concentrate on two different kinds of ar- 
rangements of the three fermions in the 3D degenerate 
Fermi gas, which configurations were also investigated in 
Ref. [14j | from the viewpoint of bipartite entanglement 
shared between two arbitrary groups of three fermions. 

(a) In the first instance a collincar arrangement is con- 
sidered, namely we put three fermions on a straight line 
numbering them in the order 1,2, and 3. The distance 
between particles 1 and 3 is r, and the intermediate par- 
ticle 2 is by a distance of x away from particle 1 (shown 
by the geometrical picture of Fig. [TJ(a)). In Fig. [TJ(a), 
the lower bound to GTE quantified by Er is plotted in 
the 3D Fermi gas according to the formula (|2Tj) in the 
function of x/r for different values kp r of the external 
fermions. The calculations can be in general performed 
only numerically, however for the limiting value kpr — * 
one obtains 



max {p tj +pjk} = P12+P23 = z -, — —, , , 2 

yfcel23,231,132 1 — x/r + (x/r) 2 



(22) 
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Substitution of this expression into (|21j) gives analyti- 
cally the curve corresponding to kpr — ► 0. Note that 
formula (|22|) holds true independently of the dimension- 
ality of the Fermi gas (i.e., both for the 2D and 3D 
cases). The curves produced in Fig. [TJ(a) exhibit two 
essential features: For any given value of fcpr, the max- 
imum of -Efj.min is achieved by the symmetrical con- 
figuration (i.e., particle 2 is located at the midpoint 
of the line connecting particles 1 and 3), still present- 
ing GTE in the system by the dimensionless distance 
kpr — 2.59. On the other, when fermion 2 starting 
from this midpoint is moved toward fermion 1 in the 
case kpr — > 0, the curve falls off to zero by the value 

x/r = 1/2(1-^3(^-2)) ~ 0.08. That is, if fcrmions 1 
and 3 are a distance r away from each other and fermion 2 
becomes closer than x ~ 0.08r to fermion 1 (or to fermion 
3 in the symmetrically equivalent situation) formula (|21[) 
does not indicate GTE among the three fermions. This 
result fits to the monogamy property of entanglement 
(28|, as in the case xjr — ► fermion 2 becomes maxi- 
mally entangled with fermion 1, excluding the existence 
of any higher order entanglement in the system. 

(b) Now let the three particles lie on the vertices of 
an isosceles triangle fermions 1 and 3 forming its base 
with length r, and fermion 2 is positioned by a distance 
of y from the midpoint of the base as it is illustrated in 
the geometrical part of Fig. [TJ(b). The curves .Ej^min in 
Fig. [2(b) are plotted against y/r for different values of 
kpr. As it can be observed, all the curves -E^mm plot- 
ted are monotonically decreasing in the function of the 
ratio y/r for any given kpr. Similarly to case (a) the 
curve corresponding to kpr — > can be treated analyti- 
cally, and one obtains vanishing GTE beyond the value 

y/r = 1/2(^3(^5-2)) ~ 0.42 (both in the 2D and 
3D Fermi gases). This supports the result of case (b) in 
Sec. Ill Dl that three fermions located at the vertices of an 
equilateral triangle (where y/r = v3/2 — 0.866) are not 
genuine tripartite entangled independent of the separa- 
tion distance r. It is also apparent from Fig.[T](b) that for 
a fixed ratio y/r, -E_R. m i n is maximal when kpr — > 0, for in 
this case the antisymmetrization effect between the three 
particular fcrmions and the rest of the fermions (which re- 
duces the amount of quantum correlations shared among 
the three fermions) becomes negligible. 



B. Fermion moving in a plane 

We now turn to the situation (pictured in the geo- 
metrical part of Fig. [2]) when fermion 2 is allowed to 
move in the two-dimensional plane given by polar coor- 
dinates (9, q) with an origin at the midpoint of the line 
of length r connecting fermion 1 and fermion 3. Let us 
restrict fermion 2 to be positioned within the circle with 
radius r/2. Recalling the definition of the GTE distance 
from Sec. U in this case the GTE distance is equal to the 
maximum separation length r between the two external 
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FIG. 1: (color online) Lower bound to the generalized robust- 
ness En of genuine tripartite entanglement shared by three 
fermions is plotted for the 3D Fermi gas for the cases (a) 
fermion 2 is moved away from the position of fermion 1 toward 
the position of fermion 3, (b) fermion 2 is moved away from 
the midpoint of the line connecting fermion 1 and fermion 3 
normal to this line (see the respective schematic geometrical 
pictures). The curves are plotted in both cases for the follow- 
ing dimensionless distances between fermion 1 and fermion 3: 
k.Fr — > (displayed in black) and fepr = 0.5, 1, 1.5, 2, 2.5, 2.59 
(displayed in color). 



fermions, below which the three fermions are genuine tri- 
partite entangled. 

In what follows, we inquire the shape of region within 
particle 2 could be located so that fermions 1,2, and 3 
by a fixed relative distance r would be genuine tripartite 
entangled. Further, owing to the symmetry of the config- 
uration it suffices to study the interval 9 £ [0,7r/2], i.e., 
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particle 2 is restricted to be situated on a quarter-disk 
of radius r/2. The polar plot of Fig. [5] represents curves 
which distinguish regions with (left side) and without 
(right side) witnessed GTE in the system, by the follow- 
ing values of the dimensionlcss distances kpr — > and 
k F r = 1,2,2.5,2.59. The witnessed GTE (i.e., to decide 
whether -ER )m m(/?3) > 0) corresponding to these curves 
for the various values of kpr was calculated according to 
the formula (|2Tj) . 

The case k F r — ► can be treated analytically yield- 
ing the result (as it can be read off from Fig. [5]) that by 
kpr — ► GTE is present in the system provided fermion 
2 is located inside a disk with radius q ~ 0.42r (gen- 
eralizing the results x ~ 0.08r and y ~ 0.42r in the 
limit kpr — > obtained in Sec. IIV Aj) . Now one can see 
from the shape of the polar curves that for greater kpr 
the corresponding disk associated with GTE squeezes to- 
ward the axis of particles 1 and 3, eventually contracting 
on the origin q = 0. It is noted that the same behavior 
would have been observed for the case of 2D Fermi g as , as 
well. This implies that if in this particular case we were 
searching for the GTE distance (i.e., the maximum sepa- 
ration r, below which GTE exists) then we should focus 
on the case where particle 2 is located at the midpoint 
between particles 1 and 3. This task will be performed 
in the sequel both for the 2D and 3D degenerate Fermi 
gases. 

V. AN UPPER AND A LOWER BOUND TO 
THE GTE DISTANCE IN THE 2D AND 3D 
DEGENERATE FERMI GASES 

A. Lower bound 

In the present section the three-fermion configuration 
is investigated, where the particles 1, 2 and 3 are posi- 
tioned evenly spaced on a straight line, with a relative 
distance r between the external particles 1 and 3. We 
develop both for the 2D and 3D Fermi gases an upper 
and a lower bound to the relative distance r, beyond 
which GTE disappears. 

Let us first calculate numerically and plot the lower 
bound to the generalized robustness of GTE defined 
by (|2"Tj) in the function of kpr. In Fig. [3] the respec- 
tive curves are plotted both for the two- and three- 
dimensional degenerate Fermi gases, and numerics shows 
that £'i?,min(/03) vanishes beyond, i.e., the lower bound to 
the GTE distance is 

r^ n = 2.3588/fc F , 

r 3 ° n = 2.5964/fc F , (23) 

respectively. Let us compare these values (|23p with the 
bipartite entanglement distance (i.e., the maximal sepa- 
ration distance between two entangled fermions), which 
are slightly smaller, and are given explicitly by the values 
[Hi 1.6163/fe F and 1.8148/fc F for the 2D and 3D Fermi 
gases, respectively. 
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FIG. 2: (color online) Witnessed genuine tripartite entangle- 
ment among three fermions is plotted for the 3D Fermi gas 
for the configuration where fermion 2 (with polar coordinates 
(8, q)) is allowed to move in a disk of radius r/2 centered at the 
midpoint between fermion 1 and fermion 3, as shown by the 
geometrical picture. The polar plot represents curves with 
black, green, cyan, magenta and brown colors, which sepa- 
rate the regions associated with GTE (left-hand side) from 
regions without having GTE (right-hand side) for values of 
the dimensionless distances k F r — > 0, k F r = 1,2,2.5,2.59, 
respectively. 

For the 3D Fermi gas the pij parameters of the state 
P3 corresponding to r^j? n are 

Pia = P23 = 0.539345, Pl3 = -0.160702 . (24) 

We can also observe in Fig. [31 as one may expect, that 
the curves are monotonically decreasing, such as in the 
bipartite case for the entanglement measure negativity 



B. Upper bound 

We continue with studying the particular three- 
fermion configuration discussed in the previous subsec- 
tion in order to establish an upper bound to the GTE 
distance, beside the lower bound already obtained. Ex- 
ploiting the mirror symmetry of the configuration, for 
any given distance r between the positions of particle 1 
and particle 3 we havepi2 = P23, thus in this case param- 
eters rfe defined by (O in Section III Al can be expressed 
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FIG. 3: (color online) Lower bound to the generalized robust- 
ness Er of genuine tripartite entanglement is plotted both for 
the two- and three-dimensional degenerate Fermi gases in the 
function of the dimensionless distance &ft", where r is the rel- 
ative distance between the positions of fermion 1 and fermion 
3, and fermion 2 is located at the midpoint, as shown by the 
upper part of the figure. 

through pi2 and px$ alone, and we obtain the relation 

r 2 = V§ri . (25) 

On the other, by applying Theorem 7 of Eggeling and 
Werner [2(| it is asserted that a three-qubit state is bisep- 
arable with respect to the partition 1 1 23 if the following 
inequalities are satisfied: 

- 1< n - 2r+ < , 
Zr 2 2 + 3r 2 + (1 - 3r+) 2 < (n - 2r+) 2 . (26) 

Now let us calculate the boundary of the area in the plane 
(n, r 2 ) described by these inequalities (f2"6")l by r + = 0.041 
and r3 = 0, which parameters correspond to the symmet- 
rical configuration of separation r^ in between fermion 1 
and fermion 3 with parameters (|24|) . The solution of the 
inequality ([2"6"]1 by r + = 0.041, r 3 = corresponds to the 
leftmost yellow shaded semi-disk in Fig. 2J representing 
states in the section r + = 0.041 and r 3 = that are 
separable with respect to the partition 1 1 23. The other 
two disks (representing biseparable states with respect to 
partitions 12|3, 13|2) can be obtained through ±27r/3 ro- 
tations around the origin of the plane {ri,r-z) [29J owing 
to the permutation symmetry of the three subsystems. 
The boundary of the convex hull of these semi-disks are 
indicated by solid blue line segments in Fig. 2J All the 
three-qubit states in the section r + = 0.041 and r 3 = 
which lie inside this polygon are in the class of bisep- 
arable states B. However, since this area corrresponds 
to a section, biseparable states in this section may exist 
outside the polygon as well. 



Next let us determine the explicit position of the 
point corresponding to the mirror-symmetrical config- 
uration with separation r^? n between the two external 
fcrmions. In this particular symmetrical configuration 
according to ([25]) the ratio r^/ri — VS and we also have 
n = (pi3 — pi2)/2. The above ratio has been displayed 
in Fig. 2] by a dashed line and the point on this line with 
coordinate ri = —0.35 (obtained by plugging the val- 
ues (l24l) into the above formula for n) corresponding to 
the separation r^? n with E R jain — has been designated 
by the red cross marker. Numerics shows, that this point 
lies outside the solid blue polygon, as it ought to owing 
to En > associated with the point. By inspection, on 
the other, this point is very close to the border of the 
polygon. 

Indeed, explicit numerical calculations yield that the 
distance r between the two outer fermions correspond- 
ing to the border of the polygon is rj^ x = 2.5988/ kp. 
This value was obtained by tuning the value of r + from 
0.041 up to ~ 0.0415 by the mean of increasing the sep- 
aration r starting from the value rj^? n so that by r^ ax 
the point represented by the the cross marker would lie 
just on the edge of the polygon. However, this distance 
r max = 2.5988/fc^ is just an upper bound to the GTE 
distance for the 3D Fermi gas. On the other, similar 
evaluations give the value r 2 ° x = 2.3599/£;f for the 2D 
Fermi gas. Comparing these values with the ones in (|23|) 
corresponding to the lower bound, it shows that the up- 
per and lower bounds to the GTE distance are indeed 
very close to each other both for the 2D and 3D Fermi 
gases. Hence, this implies that the witness of Guhne et 
al. 0] in our particular problem ought to be close to an 
optimal one. 



VI. DISCUSSION 

Previous works (e.g., (l3l. Ibj]) explored that bipartite 
entanglement exists within the order of the Fermi wave- 
length 1/k.F at zero temperature in the non- interacting 
Fermi gas and may even persist for nonzero tempera- 
tures (e.g., [lEj, l3fj). Since the system consists of non- 
interacting fermions, entanglement is purely due to parti- 
cle statistics and not to any physical interaction between 
the particles. In the present work we found the result as 
an extension of the formerly studied bipartite case that 
particle statistics is capable to generate genuine tripartite 
entanglement (GTE) as well. Furthermore, it has been 
found that the diameter of the three-fermion configura- 
tion wherein GTE is present (a lower bound to the max- 
imum diameter is given explicitly by f|23[) ) is comparable 
with the maximum relative distance between two entan- 
gled fermions, both in the 2D and 3D Fermi gases. Look- 
ing at higher order entanglement as a useful resource, 
the presence of GTE in Fermi systems would be promis- 
ing to allow for performing new quantum information 
processing tasks, exemplified by the GHZ paradox [2|. 
However, in order to do so, the amount of entanglement 
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FIG. 4: (color online). The three yellow shaded half-disks rep- 
resent areas corresponding to biseparable states with respect 
to the partitions 1|23, 12|3 and 13|2 on the (r\,ri) plane by 
the section r+ = 0.041 and r% = for the 3D Fermi gas. The 
solid blue straight lines bound the set of biseparable states 
B at this particular section. The equation ri = v3ri is rep- 
resented by the dashed line, and the red cross stands for the 
point which is located on this line with coordinate ri = —0.35. 



stored by three fermions should be somehow extracted 
from the system. In the present article, though, we did 
not consider this problem some explicit schemes has been 
put forward recently in the bipartite setting (3ll . l32l . |33| , 
some of which might be extended to the tripartite setting 
as well. 

Namely, in Ref. [3l[ it has been shown that bipartite 
entanglement can exist between non-interacting fermions 
on a lattice and can extend over multiple lattice sites even 
if the entanglement is quantified by the most restrictive 
measure, the entanglement of particles 34|. Further, con- 
sidering that in the continuum limit the entanglement of 
particles corresponds to the entanglement in the spin re- 
duced density matrix (3lT | , by continuity arguments gen- 
uine tripartite entanglement, quantified by the measure 



entanglement of particles, should exist in the lattice sys- 
tem as well. Thus, in the near future optical lattice imple- 
mentations may offer a simulation technics to observe the 
phenomenon of genuine tripartite entanglement among 
non-interacting fermions in a lattice. 

On the other, in the continuum limit the extraction 
of genuine tripartite entangled particles seems to be a 
more difficult problem: As it has been shown [l3| the 
entanglement distance between two fermions is inversely 
proportional to the Fermi momentum kp, and k F in turn 
is proportional to the density of particles. In the case of 
conduction electrons in a usual metal the density is very 
large indicating an entanglement distance of the order 
of a few angstroms. This failure might be avoided by 
using 2D electron gas formed in GaAs heterostructure, 
where the entan glem ent distance is in the order of hun- 
dred angstroms [19(] or using stored ultra-cold neutrons 
in a carefully devised experiment [33|. Also, note the 
intriguing proposal, exploiting decoherence effects to ex- 
tract bipartite entanglement created merely by particle 
statistics from semiconductor quantum wells [32|. Al- 
though the GTE distance in (j23|) is comparable (even 
greater) than the bipartite entanglement distance both 
for the 2D and 3D Fermi gases, technically these pro- 
posals appear to be very demanding when applied to the 
three-party setting. 

Finally, we would like to mention interesting future 
directions as a continuation of the present work. One 
could for example apply the same methods as in the 
present article for determining GTE distance in Fermi 
gases trapped in a harmonic trap [35| or considering GTE 
not only in spin, but in other internal degrees of freedom 
as well [36| . Also the possible existence of genuine mul- 
tipartite entanglement beyond the three-party scenario 
remains to be explored. 
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